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Abstract 

We consider the spatially homogeneous Boltzmann equation for inelastic 
hard spheres, in the framework of so-called constant normal restitution coeffi- 
cients. We prove the existence of self-similar solutions, and we give pointwise 
estimates on their tail. We also give general estimates on the tail and the 
regularity of generic solutions. In particular we prove Haff's law on the rate 
of decay of temperature, as well as the algebraic decay of singularities. The 
proofs are based on the regularity study of a rescaled problem, with the help of 
the regularity properties of the gain part of the Boltzmann collision integral, 
well-known in the elastic case, and which are extended here in the context of 
granular gases. 
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1 Introduction and main results 
1.1 The model 

We consider the asymptotic behavior of inelastic hard spheres described by the 
spatially homogeneous Boltzmann equation with a constant normal restitution co- 
efficient (see |2Ej)- More precisely, the gas is described by the probability density of 
particles f(t,v) > with velocity v G ~R N (N > 2) at time t > 0, which undergoes 
the evolution equation 

(1.1) ^ = Q(f,f) in (0,+oo)xM w , 

(1.2) /(0) = / in in R N . 



The bilinear collision operator Q(f, f) models the interaction of particles by 
means of inelastic binary collisions (preserving mass and momentum but dissipating 
kinetic energy). Denoting by e G [0, 1] the (constant) normal restitution coefficient, 



2 



when e 7^ we define the collision operator in strong formulation as 



(1.3) Q(9,f)(v)-= [ ( ^ f* ~ fg*) \u\b(u ■ a) da dv*, 

Jrn xS n-i \ e / 

where we use notations from JHj (a dual formulation shall be given in (jl.25j) which 
includes the case e = 0). Here u — v — v* denotes the relative velocity, u stands 
for u/\u\, and 'v, 'v* denotes the possible pre-collisional velocities leading to post- 
collisional velocities v,v*. They are defined by 

, v + v* 'u , v + v* 'u 

(L4) v = — + r v * = — -r 

with 'u = (1 - (3)u + (3\u\a and (3 = {e + l)/(2e) {(3 G [1, oo) since e G (0, 1]). The 
elastic case corresponds to e = 1. The function b in (jl.3|) is (up to a multiplicative 
factor) the differential collisional cross-section while B = \u\ b(u ■ a) represents the 
rate of collision of particles with pre-collisional velocities v, v * G giving rise to 
particles with post-collisional velocities G M. N defined by (jl.26j) . In the sequel 
we assume that there exists bo, b\ G (0, oo) such that 

(1.5) Vi € [-1, 1], & <&0£)< fo i, 
and that 

(1.6) b is nondecreasing and convex on (—1,1). 

Note that the "physical" cross-section for hard spheres is given by (see [TH I19p 



b(x) = est (1 — x) 



N-3 
2 



so that it fulfills hypothesis (jl.5j) and ()1.6|) when N = 3. The Boltzmann equa- 
tion (|1.1|) is complemented with an initial datum (J1.2)) which satisfies (for some 
k > 2) 

(1.7) < f in G Ll(R N ), [ f in dv = l, [ f in vdv = 0. 

JWL N JR N 

(see Subsection 11.41 for the notations of functional spaces). Notice that assuming 
the two last moment conditions in (jl.7j) is no loss of generality, since we may always 
reduce to that case by a scaling and translation argument (see |19| section 1.5] for 
instance). 

As explained in [2H], the operator (jl.3|) preserves mass and momentum: 

(L8) IL / (i) <to=0 " 

while kinetic energy is dissipated 

(1-9) j t £{f{t, ■)) = -!>(/(*,■)) 
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where the energy 8 and the dissipation functional D are given by 



£(/)=/ fWKdv, D(f):=r fU\u\ 6 dvdv*. 

Here the inelasticity coefficient r is defined by r := m 6 (^T~) anc ^ ^ e an § mar 
momentum mj is defined by 

m b := f \ b ^ . a ^ d(y= 1 r 6 ( cos gin 2 gin 7v-2 Q dQ 

JgN-l \ 2 J Jq 

(in order to get the second formula, we have set cos 8 — u ■ a). 

The study of the Cauchy theory and the cooling process of (jl.lj) - (jl.2j) was done 
in [2E| (where more general models were considered). The equation is well-posed for 
instance in L\: for < f 1VL G L\, there is a unique solution in C(R + ; L\) nL 1 (M + ; L\) 
(see again Subsection 11.41 for the notations of functional spaces) . This solution is 
defined for all times. It preserves mass, momentum and has a decreasing kinetic 
energy. The cooling process does not occur in finite time, but asymptotically in 
large time, i.e., the kinetic energy is strictly positive for all times and the solution 
satisfies 

E(t)->0 and /(£,-) 5 V=0 in M 1 (R N )-wedk * when t -> +oo, 

where M 1 (M Ar ) denotes the space of probability measures on R . We refer to 
for the proofs of these results. 



1.2 Introduction of rescaled variables 

Let us introduce some rescaled variables, in order to study more precisely the asymp- 
totic behavior of the solution. This usual rescaling can be found in [5] and ^H] for 
instance. We search for a rescaled solution g of the form 

(1.10) f{t,v)=K(t)g(T(t),V(t)v), 

where K, T, V are time scaling functions to be determined, such that K(0) = V(0) = 
1 and T(0) = (same initial datum). We choose the scaling functions K, V such 
that they are compatible with self-similar solutions, that is when g does not depend 
on time: there exists a profile function G such that 

(1.11) f(t,v)=K(t)G(V(t)v). 
In this case, the conservation of mass 



est = / f(t,v)dv = T „ N I G(w)dw 



K{t) 

v(ty 
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implies K(t) = V(t) N . The evolution equation (jl.ljl satisfied by / implies therefore 

(1.12) V'(t)V v -(vG) = Q(G,G) 

by using the following homogeneity property: for any function g such that the 
collision operator is well-defined, 

(1.13) VA e M*, Q(g(X-),g(X-))(v) = A~ (JV+1) Q(g,g)(Xv) 

(which is obtained by a homothetic change of variable). Equation ()1.12J) then implies 
that V'(t) = est =: c* > 0. When the rescaled solution g does depend on time, its 
evolution equation is 

(1.14) T\t) V(t) d t g = Q(g, g) - c* V v (vg). 

We then choose T such this equation is as simple as possible: T'(t) V(t) = 1. Hence 
we deduce the natural choice of the scaling functions 

(1.15) K(t) = (l + c*t) N , T(t) = - ln(l + ct) V(t) = (l + c*t) 

for some constant c* > 0. It is obvious that changing c* in the equation (j!.14j) only 
amounts to the multiplication of g by a positive constant and the multiplication of 
T 1 by a positive constant. In the sequel we fix without restriction c* = 1. 

Summarizing, thanks to the equation ()1.12|) and to the rescaled * variables de- 
fined by (jl.llj) . ()1.15j) . for any self-similar profile G solution to the stationary equa- 
tion 

(1.16) Q{G, G) - V« • (vG) = 

we may associated a self-similar solution F to the original equation (jl.lj) by setting 

F(t,v) = (l+t) N G((l + t)v). 
Moreover, G is obviously a stationary solution to the rescaled evolution equation 

(1-17) ^ = Q(9,9)-Vv(vg). 

which is the equation associated to ()1.1|) making the change of variables (|1.1U|) . f)l. 15|) 
(with c* = 1). Roughly speaking the re-scaling ()1 . lUj) . ()1.15|) adds an anti-drift to 
the original equation (JTTTJ) . 

More generally, for any solution g to the Boltzmann equation in self-similar 
variables 1)1. 170 . we associate a solution / to the evolution problem (jl.lj) . defining 
/ by the relation 

(1.18) f(t, v) = (l + t) N g(\n(l + t), (1 + t)v). 

Reciprocally, for any solution / to the Boltzmann equation (jl.lj) . we associate a 
solution g to the evolution problem (j!.17j) . defining g by the relation 

(1.19) g{t,v) = e- Nt f(e t -l,e- t v). 
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Given an initial datum f in = g- m e L\, we know from j2H] that there exists a 
unique solution of (jl.ljl in C(IR + , Lg) flL 1 (]R + , L\). Therefore, thanks to the changes 
of variables (|1.18|1 . (|1.19p . we deduce that there exists a unique solution g to (jl.l7j) 
in C(R+, L\) fl L 1 1 oc (lR + , L\). Moreover we have the following relations between the 
moments of / and g: 

(1.20) Vt>0, <^ 

I \\f(t,-)\-\ k \\ L i = (i+t)- k \\ g (Hi + t),-)\-\ k \\ L i. 



1.3 Motivation 

The use of Boltzmann inelastic hard spheres-like models to describe dilute, rapid 
flows of granular media started with the seminal physics paper [23I, and a huge 
physics litterature has developed in the last twenty years. The study of granular 
systems in such regime is motivated by their unexpected physical behavior (with the 
phenomena of collapse -or " cooling effect" - at the kinetic level and clustering at the 
hydrodynamical level), their use to derive hydrodynamical equations for granular 
fluids, and their applications. Granular gases are composed of macroscopic grains, 
and not microscopic molecules like in rarefied gas dynamics. The grains have only 
contact interactions, which motivates physical modelization by hard spheres with 
inelastic dissipative features built in the collision mechanism. The model of inelas- 
ticity with a constant normal restitution coefficient studied in the present paper is 
one of the simplest such model (for a more elaborated model, see for instance the 
so-called visco-elastic hard spheres model in jllj . as well as 28J). 

From the physical and mathematical viewpoint, works on the inelastic Boltzmann 
models have been first restricted to the so-called inelastic Maxwell molecules model, 
which can be interpreted as an approximation where the collision rate is replaced 
by a mean value independent on the relative velocity. Existence, uniqueness of 
solutions and the rate of decay of the kinetic energy were obtained in ["J] for the 
inelastic Maxwell molecules model with constant normal restitution coefficient. The 
Maxwell molecules model is important because of its analytic simplifications (with 
regards to the hard spheres model) allowing to use powerful Fourier transform tools. 
Polynomial tail behaviors of the self-similar profiles have been formally computed 
in |17j . Convergence to self-similarity has been established in [7J |H]. In the all, 
the inelastic Maxwell model is well understood now. Similar results have also been 
obtained for some simplified non-linear friction models in [""2*j 12*5] . 

For the inelastic hard spheres model (with constant normal restitution coeffi- 
cient), one easily sees from the discussion in Subsection 11.21 that the kinetic energy 
of self-similar solutions to equation (jl.ljl (assuming their existence) behaves like 

t^oo t z 

It is natural to expect a similar behavior for the rate of decay of the temperature 
for the generic solutions to equation (jl.ljl . This conjecture was made twenty years 
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ago in the pioneering paper [21], and this rate of decay for the temperature is 
therefore known as Haff's law. Such a law is a typical physical feature of inelastic 
hard spheres which does not hold for inelastic Maxwell molecules model (for which 
the temperature follows an exponential law). Let us emphasize that in jB] a pseudo- 
Maxwell molecules model was considered (multiplying the collision operator by some 
well-chosen scalar function of time), restoring Haff's law and still preserving the nice 
simplications of Maxwell molecules. However some important aspects of inelastic 
hard spheres such as the tail behavior are not preserved by this model. 

On the basis of the study of the particular case of Maxwell molecules, Ernst and 
Brito JH] conjectured that self-similar solutions, when they exist, should attract any 
solution, in the sense of convergence of the rescaled solution. They also conjectured 
and formally computed some "over-populated" tail behaviors for the self-similar 
profile (namely decreasing slower than the Maxwellian) depending on the collision 
rate (see [T7j for instance). 

Recently the works ^IJJ E] laid the first steps for a mathematical analysis of the 
inelastic hard spheres model (with constant normal restitution coefficient). In [TU] . 
the authors proved the existence of steady states and gave estimates showing the 
presence of over-populated tails for diffusively excited inelastic hard spheres, that is 
when one provides an input of kinetic energy to the system preventing the collapse, 
which is modelized by some A„/ term added to equation (jl.lj) . In jU] a priori integral 
estimates on the tail of the steady state (assuming its existence) were established 
for the spatially homogeneous inelastic Boltzmann equation with various additional 
terms, such as a diffusion, or an anti-drift as in (|1.17|) . 

In the present paper, we prove, for spatially homogeneous inelastic hard spheres 
with constant normal restitution coefficient, the existence of smooth self-similar 
solutions, and we improve the estimates on their tails of jU] into pontwise ones. We 
also give a complete regularity study of the generic solutions in the rescaled variables, 
as well as estimates on their tails. In particular, we give the first mathematical proof 
of Haff's law and we show the algebraic decay of singularities. 

In a forthcoming work p9, , we shall prove the uniqueness and the local stability 
of these self-similar solutions for a small inelasticity. 



1.4 Notation 



Throughout the paper we shall use the notation (•) = ^1 + 
any gGl, the Banach space 



2 . We denote, for 




More generally we define the weighted Lebesgue space 
q G M.) by the norm 





i i/p 
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when p < +00 and 

||/||l«(r*) = supess^ 6MW (v) q 

when p = +00 (where supess denotes the essential supremum). 

The weighted Sobolev space W£' P (M. N ) (p G [1, +00], g G R and fc G N) is defined 
by the norm 



1/p 



|s|<fc 



where (9 s denotes the partial derivative associated with the multi- index s G N^. In 
the particular case p = 2 we denote if* = W 7 ^' 2 . Moreover this definition can be 
extended to if* for any s > by using the Fourier transform. 

We also denote by L} oc (£l) the space of locally integrable functions on a given 
set C M. N , that is the space of measurable functions on f2 which are integrable on 
every compact subset of Q. Finally for h G M. N , we define the translation operator 
T h by 

Vv G R N , T h f(v) = f(v-h), 

and we shall denote by "C" various constants which do not depend on the collision 
kernel B. 

1.5 Main results 

In this subsection we consider a normal restitution coefficient e G (0, 1) (except in 
dimension N = 3 where the case e = can be included, see the discussions in the 
proofs) . 

First we state a result of existence of self-similar solutions. 

Theorem 1.1 For any mass p > 0, there exists a self-similar profile G with mass 
p and momentum 0: 

0<GeLl Q(G,G) = V v -(vG), J G ( M dv=(?Y 

which moreover can be built in such a way that G is radially symmetric, G G C°° 
and 

Vv G R N , aie- a2H < G(v) < A^'^ 
for some explicit constants ai,a,2,Ai,A 2 > 0. 

Second we establish that Haff 's law holds. 
Theorem 1.2 For any p G (1, +oo) ; p > 0, and some initial datum f in such that 



0</inG^nF, / f-J) 



dv 
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the associated solution of the Boltzmann equation M.lM.ty) in C(R + ; L^flL 1 (R + ; L\) 
satisfies Haff's law in the sense: 

(!-21) Vt>0, rz<£(t)< 



;i + t) 2 - w - (i + t) 2 

for some explicit constants m, M > depending on the collision kernel and the mass, 
kinetic energy and LP norm of f in . 

Third we give a more precise and general result on the regularity and asymptotic 
behavior of the solutions g to the rescaled equation (jl.l7|) (note that the first point 
of Theorem 11.31 implies in particular the preceding theorem). 

Theorem 1.3 For any p G (1, oo), p > 0, and some initial datum g m such that 



0<g- m eLlnlS, [ 9in(l) 
Jrk V v J 



dv 



the unique solution g in C(R + ]L\) fl L 1 1 oc (M + ; L\) of \l.l r l\j with initial datum g- m 
satisfies: 

(i) It remains bounded in LP for all times, with uniform bound as t goes to infinity. 
Similarly for any q > 0, if gi n G LP q , then the solution remains bounded in LP q 
for all times with uniform bound as t goes to infinity. As for the Sobolev 
norms, for any s,q > 0, there is w > such that if g m G H q+W , then the 
solution remains bounded in iff for all times, with uniform bound as t goes to 
infinity. 

(ii) For any arbitrarily large s,q>0, there exists A > and some decomposition 
g = g s + g R of the solution g such that g$ > and 

sup ||<?f || < +oo, \\g*\\ j = O {e- xt ) . 

(Hi) Concerning the tail behavior, we have the following lower and upper bounds: 
there are some explicit constants a 1; a2 > such that 

VveR N , \immig(t,v)>a 1 e~ a ' M , 

and, for any r > and s G [0, 1/2), there are some explicit constants Ai,A 2 > 
such that ( appearance of exponential moments ) 



Vt>r, [ g(t,v)e~ AlMs dv < A 2 . 



Moreover similar integral upper bounds with s G [1/2, 1] are uniformly propa- 
gated in time if they are satisfied for the initial datum. 
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All the constants in this theorem can be computed in terms of the mass, kinetic 
energy and the different norms assumed on g in , and the parameters. 

Remark 1.4 Note that point (ii) of Theorem implies, by coming back to the 
original variables, that for an initial datum < f\ n G L\ R LP, p G (l,oo), the 
unique associated solution of M.l)) in C{R + ]L\) PI L 1 (M. + ;Ll) satisfies a similar 
decomposition as above, but where the remaining part decreases with polynomial and 
not exponential rate. Hence we have shown that the amplitude of the singularities 
decreases algebraically in the original variables. 

More precisely it is likely that, in the original variables, singularities far from 
decrease in fact exponentially fast, whereas those close to cannot decrease faster 
than polynomially (due to the fact that the damping effect of the loss part of the 
collision operator degenerates at this point). 

1.6 Method of proof 

The main tool in this paper is the regularity theory of the collision operator: we 
show that its gain part satisfies similar regularity properties as in the elastic case 

ESI EDI 123 E2] • Following the study in the elastic case in (32j , we deduce uniform 
propagation of Lebesgue norms for the solutions in the rescaled variables ()1.17j) . 

A first consequence is that the temperature in the rescaled variables is uniformly 
bounded from below by some positive number as soon as the initial datum satisfies 
some LP bound. Translating this estimate in the original variables, it proves Haff 's 
law. 

A second consequence is the existence of self-similar profiles (or steady states) 
for the rescaled equation (|1.17|) . which provides existence of self-similar solutions 
for the original equation This existence result is proved by the use of a 

consequence of Tykhonov's fixed point Theorem (see Theorem 14. which is an 
infinite dimensional (rough) version of Poincare-Bendixon Theorem on dynamical 
systems, see for instance [31 Theoreme 7.4] or [THIHE]- ft states that a semi-group 
on a Banach space y with suitable continuity properties, and which stabilizes a 
nonempty convex weakly compact subset, has a steady state inside this subset. We 
apply it to the evolution semi-group of (|1.17|) in the Banach space y = L\. The 
existence and continuity properties of the semi-group were proved in |2E] and the 
nonempty convex weakly compact subset of nonnegative functions with fixed mass 
and momentum and bounded moments and LP norm, p G (1, +oo) (for some bound 
big enough) is stable along the flow thanks to the above uniform LP bounds in the 
rescaled variables. 

Still following the study in the elastic case in [22], we also deduce from the 
regularity properties of the collision operator the uniform propagation of Sobolev 
norms as well as the exponential decay of (the amplitude of) the singularities for 
the solutions in the rescaled variables ()1.17|) . That straightforwardly implies the 
smoothness of self-similar profiles as well as the algebraic decay of (the amplitude 
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of) singularities for solutions to the Cauchy problem in the original variables (jl.lj) - 

Let us now turn to the study of the tail behavior. On the one hand, we prove 
lower pointwise estimates on the self-similar profiles by the mean of some elemen- 
tary maximum principles arguments inspired from |19j . On the other hand, we 
prove explicit lower pointwise estimates on generic solutions in self-similar variables 
using the spreading effect of the evolution semi-group associated to (|1.17j) (in the 
spirit of [T3J EH1 EI])- Finally, upper pointwise estimates on the self-similar pro- 
files are obtained using moments estimates established in [5| and elementary o.d.e. 
arguments. 



1.7 Weak and strong forms of the collision operator 

Under our assumptions on b, the function a 1— > b(u ■ a) is integrable on the sphere 
S^ -1 , and we can set without restriction 

b(u ■ a) da = \§ N - 2 \ / b(cos9) sm N - 2 9d9 = l. 
-1 Jo 

Thus we can write the classical splitting Q = Q + — Q~ between gain part and loss 
part. The loss part Q~ is 

(1.22) Q-(g,f)(v):= ( [ g(v*)\v-v m \dvA f(v) = (g*$)f, 

\Jr n / 

where $ denotes $(z) = \z\. For any distribution g satifying the moment conditions 
J RN g dv = 1, f RN g v dv = 0, we have (see for instance j2Hl Lemma 2.2]) 

(1.23) (#*$)> |<4 
The gain part Q + is defined by 

(1.24) Q+(g,f)(v):= [ 11^1 \ u \b(u ■ a) da dv*. 

In the sequel, we shall need two other representations. On the one hand from 
[Hj, there holds: for any ip £ Lf, f,g £ L\ 

(1.25) / Q + (g, f)(v)t/j(v) dv = / / g* \u\ b(u ■ a) ip(v') da dv* dv, 

where v' denotes the post-collisional velocity defined by 

/ v + v* u' . 1 — e 1 + e. . 

(1.26) v = H , u = u-\ \u\a. 

K ' 2 2' 2 2 1 1 

On the other hand, we shall establish a Carleman type representation for granular 

gases: 
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Proposition 1.5 Let E%, v be the hyperplan orthogonal to the vector v — 'v and 
passing through the point Q(v,'v), defined by 

Q(v, 'v) := v + (1 - /T 1 ) (v - 'v) = (2 - /T 1 ) v + (/T 1 - 1) 'v. 

Then we have the following representation of the gain term 

nN-l r r 

(1.27) Q+( g J)( v ) = —— / / \ v - v ^-"B\'v-v\- x 'g*'fd'vdE('v # ). 



aN-\ p 2 I / 1 
P e J'vm N J'v,eE e , 

v , '1 

Recall that (3 = (1 + e)/ (2e) and B := B(u, a) = \u\ b(u ■ a) . 

Remark 1.6 Let us emphasize that in dimension N = 3 the expression \1.21\l of 
the gain term Q + makes sense even when e = (since (e(3) 2 converges to 1/4 and 
Q(v, f v) = 2v — 'v for e = 0), while the formula defining Q + in M.ty) seems to be 
singular when e — > 0. Hence this Carleman representation allows to define a strong 
formulation of Q + for e = 0, at least in the physical case of the dimension N = 3. 

Proof of Proposition 11.51 We start from the basic identity 

(1.28) - / F(\u\a - u) da = , * / 5(2 x ■ u + |x| 2 ) Fix) dx, 
2 Jsiv-i \u\ n J r n 

which can be verified easily by completing the square in the Dirac function, taking 
the spherical coordinate x + u = r a and performing the change of variable r 2 = s. 
We have the following relations from (jl.4j) 



' v = v + (P/2) (\u\a-u) 
'u* = - ((3/2) (\u\a - u) 



(1.29) 

and thus starting from the strong form of Q + we get 

Q+(g,f) = e- 2 [ Bf(v + ((3/2) (\u\a - u) ) gU-((3/2) (\u\a - u)) dv* da. 

Applying (jl.28|) yields 

Q+(gJ) = 2e- 2 [ \u\ 2 ~ N B5(2x-u+\x\ 2 ) f(v + (f3/2)x)g(v*-((3/2)x)dv Jf dx. 

JR N xR N 

We do the change of variable x — > 'v = v + ((3/2)x (with jacobian ((3/2) N ). Then, 
keeping 'v fixed, we make the change of variable v* — > 'v* (with jacobian 1 since 
'v * = v + v * — 'v) . This gives 

oN+i r 

1 1 |2 — JV D r/o ™ „. i 1™^ 



Q + (9, f) = ^n \u\ 2 ~ N BS(2x ■ u + |x| 2 ) f{'v)g{'v.) d'v*d'v. 

P e JR N xR N 
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Finally, keeping 'v fixed, we decompose orthogonally the variable 'v* as v + V% n + V<i 
with V\ — ('f* — v) ■ n, n — ('v — v)/\'v — v\ and V 2 orthogonal to {'v — v). Let us 
compute the Dirac function in the new coordinates. Since x = (2/(3) ('v — v) and 
u — (v — f*), 

2x-u+\x\ 2 = (4/f3)('v-v)-(v-v*) + (4/f3 2 )\'v-v\ 2 

= [A/ (3) ((/T 1 -l)\'v-v\ 2 + ('v-v)-('v-vS 

From the momentum conservation ('v — v*) = (v — 'v*) and the orthogonal decom- 
position above: 

2 x ■ u + \xV 



(4/ 73) - 1) | 'u - *;| 2 - K | 'v - u|) . 



Hence we obtain the following representation: 

2 N+1 



Q + (g,f) 
5 



(3 N e 2 
A\'v - v\ 



\u\ 



\2-N 



B 







xR JV - 1 xK JV 

{(3~ x -l)\v-v\ 



V-i 



f('v)g(v + Vm + V 2 ) dV x dV 2 d'v. 



It remains to remove the Dirac mass: we use the obvious identity 



5 



A\'v - v\ 



(/T- l)\'v -vl-Vt F{V x )dVj 



(3 



4 \'v — v 



F({(3- x -l)\'v- 



to finally obtain representation (jl.27|) . □ 

The parametrization by the Carleman representation means that for v and 'v 
fixed, the point 'v* describes the hyperplan orthogonal to ('v— v) and passing through 
the point Q(v, 'v) on the line determined by v and 'v. Note that in the elastic case, 
Q(v, 'v) = v, whereas here Q(v, 'v) is outside the segment [v, 'v], which reflects 
the fact that for the pre-collisional velocities, the modulus of the relative velocity is 
bigger than \v — v*\. In the limit case e = 0, Q(v, 'v) = 2v — 'v. 

The geometrical picture (in a plane section) is summerized in Figure ^ 
From this proposition we immediately deduce the following representation, which 
is closer to the classical Carleman representation for the elastic Boltzmann collision 
operator. From (jl.29j) we deduce 

J3_ 
V2 



v — v\ 



\u 



y/l-(u-a). 



Hence we get 

Q + (f,g)(v) 

with 

a 



a 



2^ 



7 



l v \N-2 



'g* b(u ■ a)d'v* > d'v, 



and 



b{x) 



x 



-(N-3) 



b(x). 
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n 



V 

-*- 



V* 



Figure 1: Carleman representation for granular gases 



2 Regularity properties of the collision operator 

In this section the final goal is to estimate quantities such as 

for p > 1, i.e., the action of the collision operator on the evolution of the LP norm 
(to the power p) of the solution along the flow. We shall use minoration estimates on 
Q~ deduced from (jl.22|) - (|1.23|) . together with convolution and regularity estimates 
on Q + . The latters seem to be new in the inelastic framework but they are an 
extension of similar estimates in the elastic case e = 1. 

The estimates on Q + can be splitted into tree groups. First the convolution-like 
estimates, which originated (in the elastic case) in the works of Gustafsson (2H 122] 
(see also [221113123) and were first extended to the inelastic case in jTHj. Second 
the regularization estimates in Sobolev spaces which originated in the elastic case 
in the works of Lions [2H], Bouchut and Desvillettes [TU], Lu (211 ( see a l so E2] for 
some extensions). Third the non-concentration estimates in L , which originated in 
the elastic case in the work of Mischler and Wennberg [30J, and were extended by 
Abrahamsson [I]. Let us also mention that regularity properties of Q + are reminis- 
cent of the work of Grad on the linearized collision operator |23j . The main tool 
to extend the second group of estimates to the inelastic case shall be the Carleman 
representation for granular gases of Proposition 11.51 (the third group of estimates 
can be extended with this tool as well, see [29J). Before turning to the regularity 
study of Q + , we recall convolution-like estimates. 
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2.1 Convolution-like estimates 

In the elastic case e = 1, convolution-like estimates for the gain part of the collision 
operator were first proved in [^TJ|22]. This proof was simplified by a duality argu- 
ment in [32], where also a more precise statement was given. These estimates were 
extended to the inelastic case, for a constant normal restitution coefficient e G [0, 1], 
in (i n a form slightly less precise than in [22]). Also a result weaker in one 
aspect (less precise for the treatment of the algebraic weight) but more general in 
another (valid in any Orlicz spaces, and valid for more general collision kernels) was 
proved in [2E] ■ Here we only state the precise result we shall need, whose proof is 
straightforward from the arguments in jTT^ Proof of Lemma 4.1] and (S3 Proof of 
Theorem 2.1]. 

We make the following assumption on the cross-section: no frontal collision 
should occur, i.e., b(cos9) should vanish for 9 close to n: 

(2.1) 39 b >0; support b (cos 9) C {9 / < 9 < tt - 6 b } . 

To exchange the roles of / and g, we introduce the symmetric assumption that no 
grazing collision should occur, i.e., 

(2.2) 39 b >0; support b (cos 9) C {9 / 9 b < 9 < vr} . 

Then we have (from the proofs of ^01 Lemma 4.1 and Proposition 4.2]): 

Theorem 2.1 Let k,t] G M., p G [1, +oo], and let B = $6 be a collision kernel with 
b satisfying the assumption \2.1\) . Then for any e G [0, 1] the associated gain term 
satisfies the estimates 

\\Q + {gJ)\\ LP <C k ^ p {B)\\g\\ L] M ||/|U , 

11 11 \k+7]\ + \r]\ k + ri 

with 

C k , v , p (B) = C (sm(9 b /2)r^)-yp' ||fo|| il(s „_ X) ||$|| L . . 

If on the other hand assumption \2.1\) is replaced by assumption \2.2}) . then the same 
estimates hold with Q + (g, /) replaced by Q + (f,g)- 

2.2 Lions Theorem for Q + 

In this subsection we assume that the collision kernel B = $ 6 satisfies 

(2.3) $ G C^°(R JV \ {0}), fcGCn-M). 
Then we have the 

Theorem 2.2 Let B be a collision kernel satisfying \2.'J\i . Then for any e G (0, 1], 

the associated gain term Q + satisfies for all s G K.+ and r\ G M+ 

\\Q + (gJ)\\ H s + (N-^<C(s,B) IMI^II/II^ 

for some explicit constants C(s, B) > depending only on s and the collision kernel. 
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Remark 2.3 In dimension N = 3 this theorem extends to the case e = 0, with 
uniform bound C(s, B) for e G [0, 1]. The only obstacle to the treatment of the case 
e = indeed is the constant 2 N ~ l (3~( N ~ l ) e~ 2 in front of the Carleman representa- 
tion, which may blow up as e —> 0. 



Proof of Proposition 130 We follow closely the proof of inspired from the works 
of Lions and Wennberg [35] . Indeed the Carleman representation proved above 
in Proposition 11.51 allows essentially to reduce to the study of the elastic case. 
We assume first that r\ = 0. We denote 

B(\'v- 'v,|,|'t;-t;|) 



B(\v- 


- 


, COS 9) 


\v — u*| 


N-2 


| 'v — v\ 



which belongs to C^°((1R + \ {0}) 2 ) under assumption (|2.3|) . We define the following 
(Radon transform type) functional: for g smooth enough, Tg is defined by 



Tg{y) = / B(z,y)g(z)dz 

J ny+y x 

with /i = (2 — Let us relate this functionnal with the Carleman representa- 

tion (jl.27j) : we have 



and since 
we deduce 



/ I 1 2— N d 1/ 1-1 ' j / 

/ \v — v*\ B | v — v\ g*d v* 

J'v,£E e , 

B{\'v*—'v\,\v—'v\) 'g*d'v* 

'«» e Q(v, 'v)+(v- 'v) 1 - 

B(\z\, \v - 'v\) r_ /^(z) dz, 



'v- t £E e , 



v - v*\ 2 N B\'v — v\ 1 'g*d'v* 
B(\z\, \v - 'v\) r_ 

fj-(v— 'v)+(v— 'v) 1 - 

(t, v oTor_ ,„) (^)(v). 



Hence the representation (jl.27j) writes 



2 AT-I 

/3 1 



G + fo>/)(») = 7ijv=T72 / f{'v){T, v oTo T - lv ){g){v)d'v. 
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Thus if one has a bound on T of the form 

(2.4) \\Tg\\ Hs+{ N-w<C T \\g\\ H *, C T > 0, 

then by using Fubini's and Jensen's theorems one gets 

\\Q + (gJ)\\ 2 HsH N-W < C\\f\W [ f('v) \\(T, v oToT_, v )(g)\\ 2 HS+iN _ 1)/2 d'v 

< C\\f\\ L i [ fCv)\\(ToT- lv )(g)\\ 2 HS+{N - 1)/2 d'v 

Jm N 

< C7C r ||/|| L1 / f('v) \\T-, v g\\%, d'v 

< CC T \\g\\ 2 HS \\f\\# / f('v)d'v<CC T \\g\\ 2 Hs \\f\\ 2 L1 , 

which concludes the proof. Thus it remains to prove (|2.4|) . But, up to an homothetic 
factor, T is exactly the operator which was studied in detail in |J5] and [32]. More 
precisely, 

Tg(y) = fg(/iy) 

where T is the Radon transform 



Tg(y) = / B(z,y)g(z) dz, 

Jy+y 1 - 

which was introduced in the elastic case in associated with a kernel B related 
to our collision kernel by 

B(z,y) = B(z, f jL- 1 y). 
It was proved in [^21 Proof of Theorem 3.1] that 

||TV/||ir»+c*r-i)/a < C \\g\\ H . 

for an explicit bound C depending on some weighted Sobolev norms on B. Coming 
back to T, we obtain ([2.4)1 . This ends the proof when rj = 0. The extension to rj > 
is straightforward (and exactly similar to |3*2l Proof of Theorem 3.1]). □ 

As a Corollary we deduce from Theorem 12 . 21 the following estimate in Lebesgue 
spaces by Sobolev embeddings (the proof is exactly similar to p2, Proof of Corol- 
lary 3.2]). 

Corollary 2.4 Let B be a collision kernel satisfying \2.§i . Then, for all p G 
(1, +oo) ; rj G K ; we have 

\\Q + (gJ)\\ LP <C( P , v ,B) \\g\\ L ,\\f\U t 

11 11 'I 2|?7| 

where the constant C(p, rj, B) > only depends on the collision kernel, p and r\, and 
q < p is given by 

{2N ~ 1)P if P e(l,2N) 



(2.5) 



N+(N-l)p 
^ tfpe [2iV,+oo) 
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2.3 Bouchut-Desvillettes-Lu Theorem on Q + 

Now we turn to a slightly different regularity estimate on Q + , which is a straightfor- 
ward extension of the works fOJ|2I] in the elastic case e = 1. This class of estimate 
is weaker than Lions's Theorem 12.21 since the Sobolev norm of Q + is controlled by 
the square of the Sobolev norm of the solution with smaller order, which does not 
allow to take advantage of the L 1 theory. Nevertheless, it is more convenient in 
other aspects since it deals directly with the physical collision kernel. 

Theorem 2.5 Under the assumptions made on B in Subsection for any e G 
[0, 1] the associated gain term Q + satisfies, for all s G M+ and rj G M. + , 



\\Q + (g,f)\\ H s«N-iy2<C(s,B) 



\9\\H, 



17+2 



H 



V+2 



\Q Wr 1 



r 1 

%+2 



for some explicit constant C(s, B) > depending only on s and B. 



Proof of Theorem \2. 51 We follow closely the method in [TP] . We write it for rj = 
but the general case is strictly similar. 

Let us denote F(v, u*) = f{y) g(v*) \v — v*\. The same arguments as in [TU] easily 



lead to 



where TQ + denotes the Fourier transform of Q + according to v, F denotes the 
Fourier transform of F according to v,v*, and 



c+ 3 ~ e c , 1 + e Id 



Thus 



\FQ + (0\ 2 <\\b\\l H s»-i) 
Let us consider frequencies £ such that |£| > 1. As 

\ne,o\ 2 da 



\F(e,C)\ 2 da 



S JV-1 



S^- 1 J\z\ 



+oo 



d_ 

dr 



3-e 1 + e 1 + e 

c H re, c 

4 4 4 



1 + e 



■ ro 



da dr 



< C 



+oo 



l£l 



r , 3-e 1+e 1+e 1+e 
r I — : — q H — ra, £ : — ro ~ 



< C 



(V 2 - Vi)F 
F 



4 

3-e 
4 



4 



4 



.1+e 1+e 1+e 
£ H ^ — rcr ) — ^ — £ : — rcr 



4 



3-e l + e r 1 + e 1 + e 



IC|>KI 
(V 2 -Vi)F 



4 



4 



3-e 1+e 1+e _ 1+e 



4 



4 



4 



4 



ICI 



JV-1 ' 



18 



where we have made the spherical change of variable £ = rcr, we deduce 



+ /t\|2 |t|2s+(AT-l) 



ICI>i 



i<iei<lci 



Fl iz£ £ + I+£ ( ,I+£ £ _I+£ C 



4 



4 



x 



/ 3 — e . 1 + e^l + e^, 1 + e A 
(V 2 - Vi)F( — ^C— 



4 - 4 " 4 
Finally we make the change of variable 



|t|2a+(JV-l) 



4 * 4 " 4 ' 4 

(whose Jacobian is uniformly bounded from above and below for e G [0, 1]) to obtain 



\^Q + m 2 \e s+{N - 1] ^ < c^ 2 



l€l>i 



F(X,F) 



x 



(V 2 - Vi)F(X,Y) (X} 2s {Y) ls dX dY 



< C Wbf^-^W u- ||( r-r : )t || „- 
12 



< CUftll^-xjIbll^ll/lli.. 

Then small frequencies are controlled thanks to the L 1 norms of / and g, which 
concludes the proof. □ 



2.4 Estimates on the global collision operator in Lebesgue 
spaces 

We consider a collision kernel B = $6 with $(«) = \u\ and b integrable. We shall 
make a splitting of Q + as in [321 Section 3.1]. We denote by 1e the usual indicator 
function of the set E. 

Let : K — > R + be an even C°° function such that support G C (—1,1), and 
j R Qdx = 1. Let : — > M + be a radial C°° function such that support© C 
5(0, 1) and J RN Qdx — 1. Introduce the regularizing sequences 

{0 m (i) = mO(mx), iGl, 
e n (ar)n iV e(na:), a; G 

We use these mollifiers to split the collision kernel into a smooth and a non-smooth 
part. As a convention, we shall use subscripts S for "smooth" and R for "remainder" . 
First, we set 

$S,n = ©n * ($ IaJ , = * ~ $S,n, 
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where A n stands for the annulus A n = {16 M. N ; ^ < \x\ < n}. Similarly, we set 



bR,m = 6-6 



S.mi 



where Z m stands for the interval Z„ 



{s G 



< Id < 1 



is un- 



derstood as a function defined on R with compact support in [—1,1]). Finally, we 

set 

Q = Q s Qri 

where 

®s,n ( I v - I ) 6 5 ,m (cos 0) '/ da dv* 



Qt(gJ) = e~ 2 [ 



K JV x giV-l 



and 



Qr — Qrs + Qsr + Qrr 



with the obvious notation 

QR S (gJ) = e- 2 



$R,n bs, m 'g* 7 dv* da 



R N xS N-l 



QUg,f) = e' 2 [ 

JR 

QU9,f) = e- 2 [ 



T xS N ~ 1 



®s,n b R , m 'g* 7 du* da 



®R,n b R , m 'g* 7 d«* do"- 



r jv x §jv-i 

Now we follow the proof as in [H2J Section 4.1] since we have the same functional 
inequalities in Sobolev and Lebesgue spaces, using also some ideas from [T^l [2H] to 
simplify it. 

Proposition 2.6 Let us consider e G (0, 1] and the associated gain term Q + . For 
any e > 0, there exists 6 G (0, 1), only depending on N and p, and a constant 
C e > 0, only depending on N, p, B and e (and blowing up as e — > 0), such that 



[ Q + (fJ)f P ~ 1 dv<C £ 



l+p9 n ^||p(l-0j 
LP 



1/p 



Remark 2.7 For general dimension N > 2, the estimates in this section are valid 
only for e G (0, 1] (the constants may blow up as e 0). However when N = 3 
they are uniform on e G (0,1] and then extends to the limit case e = (see also 
B,emark \1.6]) . 



Proof of Proposition VJM Let us fix e > 0. We split Q + as Q$ + Q RS + Qs R + Q RR 
and we estimate each term separately. Remember that the truncation parameters n 
(for the kinetic part) and m (for the angular part) are implicit in the decomposition 
of Q+. 
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By Corollary 12 .41 there exists a constant C(m, n) > 0, blowing up as m or n goes 
to infinity, such that 

\\QtU\f)\\ LP <C(m,n) U/M/IUi, 
for some q < p defined in f!2.5|) . Hence by Holder's inequality, 



F- l Qt{fJ)dv< 



Fdv 



p-i 



(Q+ydv 



< 



p-i 

LP 



\\Qt(f,f)\\ LP <C(m,n) 



L"\\.I\W II/Hlp 1 - 



Next we fix a weight rj > — 1 and we estimate the norm of Qsn(f,f) and 
Qimifif)- We use that ||2>j2 jm ||x, 1 (s JV '— x ) S oes to as m goes to infinity (since b is 
integrable on the sphere), and we obtain, using Theorem l2.1l with k = 1 and splitting 
the angular integration between a part u ■ a < with no grazing collision and a part 
u ■ a > with no frontal collision, 



WQUfJ), Q%R(f,f)\y<<rn) 



1 1 + r/ 1 + | -r] | 



for some e(m) going to as m goes to infinity. Since 1 + 77 > 0, we can write 
|1 + 77 1 + \rj\ — 1 + 2?7 + , where 77+ = max(r7, 0). Hence by Holder's inequality, 



/ f p - 1 (Qt R (f,f) + Qt R (f,f))dv= I (f(v)^ 

Jm N v J Jm N 



l/p\p-l ^CSR 



1 QsR + Q 



RR 



< 



< 
p-i 



(Q + SR + Q + RR)(v)- W ) P dv 



% [\\QUfJ)\\L* 1/t + \\QUfJ)\\L* 1/t )<e 

l/p V -!/P -1/P' / 



l/p 



It remains to estimate the term corresponding to Qr S - We have the trivial 
estimate 

®R, n <Cn~ l (\v\ 2 + \v*\ 2 ) 

from which we deduce that 



Jm N 



< On' 1 / b s , m f U(fy- l (\v\ 2 + \v*\ 2 )dvdv*da:=h + I 2 . 

JR N xR N xS N 

Now we treat separately the two terms of the right-hand side: 



h 



bs,m{f\v\ 2 )U{f) p - l dvdv* da 



R N xR^xS" 

< C [ b s ,m IF + {f') p ] (f \v\ 2 ) dvdv* da := I M + I 1>2 

JR N xR N xS N 
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[2.7) J 1>2 = C I I b s , m KJ J - f\v\ 2 dv'dvda<C(m)\\f\\ Ll iij UP . 

b*(0) J ° (p [V ) 



using Young's inequality xy p 1 < (1/p) x p + ((p — l)/p) y p on the product /*(/') p 1 - 
The control of Ixi is immediate by integrating separately the angular variable: 

/i,i<c||/IUi II/IIlp. 

For Ii t 2, as in the proof fl9\ Proposition 4.3] and using the notations of |2EJ Lemma 
4.4], we make the change of variable i>* i— > v ' — 0*(t>*) = 4>l, v ,a( v *) keeping v, a 
fixed, which is a C^-diffeomorphism from O = {t>* G M. , ra/l} onto its image. 
Thanks to [2E1 Lemma 4.4] and because bs, m has compact support in (—1,1), its 
Jacobian J* = det(D<p*) satisfies 

(2.6) C{mY l < J*(>*) < C(m) Vw* G u ■ a £ support 6 S>TO , 

for some constant C(m) G (0, oo) which blows up when m goes to oo. Hence we 
straightforwardly deduce 

iff 

'S,m ~fTZ 

xS N J<f>*{0) 

The term J 2 is treated in a similar way: it is splitted as above using Young's in- 
equality on the product f(f') 1P ~ l - The term J 2 i involving f p is directly estimated as 
for the term For the term J 22 involving (f') p we proceed as for the term 7 12 . 
We make now the change of variable v 1— > v' — <f> e ,v,a( v ) keeping v*,cr fixed, where 
we use again the notations and results of [2E1 Lemma 4.4]. Since its Jacobian J 
also satisfies the bound (|2.6|) . we get the estimate ()2.7|) for the term J 2)2 . We finally 
deduce 

/<^ll/IUj ll/IU- 

n 2 

Gathering the previous estimates we deduce 

F-'Q+iLfidv < C(m,n) H/llMll/ll^ll/lli; 1 

+^ll/llx* \\fh P + e(m) \\f\\ Ll \\f\\l P 
n 2 2 ^i/p 

where q is defined by (12.5)1 . e(m) goes to as m goes to infinity, and C(m, n), C(m) > 
0. Hence for any given e > 0, by first fixing m big enough, then n big enough, we 
get 

/ f p - l Q + U, f) dv < C £ \\f\\ Lq \\f\\» \\f\\% 1 + e \\f\\ Ll H/IIJ, 

for some explicit constant C £ > 0. Combining this with elementary interpolation, 
we deduce that there exists 9 G (0, 1), only depending on N and p, and a constant 
C £ > 0, only depending on N, p, B and e, such that 

/ f p - l Q + UJ)dv <C7 e ||/||^||/||i^||/||£: l + e ||/mi|/||^ 

Jr n ' x lv 

1/p 



This concludes the proof. □ 
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3 Regularity study in the rescaled variables 



In this section we show the uniform propagation of Lebesgue and Sobolev norms 
and the exponential decay of singularities for the solutions of (jl.l7|) . 



3.1 Uniform propagation of moments - Povzner Lemma 

Let us prove that the kinetic energy of g remains uniformly bounded from above as 
t goes to infinity. Using (jl.l7J) and (jl.9|) . we get 



— / g\v\ 2 dv < —t I g g* \u\ 3 dv*dv + 2 / g\v\ 2 dv. 

R N xl" 



dt 



On the one hand, from Jensen's inequality (see for instance [23 Lemma 2.2]), there 
holds 

g* \u\ 3 dv* > p \v\ 3 . 



On the other hand, Holder's inequality yields 

1/3 / f \ 2/3 



/ g\v\ 2 dv< I gdv\ I g\v\ 3 dv I , 
Jr n \Jr n J \Jr n J 



which implies that 

/ g \v\ 3 dv > p 1/2 ( [ g\v\ 2 dv\ . 
Jr n \Jr n J 

Thus 

^- ! g \v\ 2 dv < -t p 3/2 ( [ g\v\ 2 dv) + 2 ( [ g\v\ 2 dv 
dt Jrn \Jrn J \Jrn 

< rp 3/2 I / g\v\ 2 dv 



T p3/2 



- 1/2 

9 \ v \' 2 dv 



and by maximum principle we deduce 



(3.1) sup / g \v\ 2 dv < Cjjmax <j ( — ^- ) , / g in \v\' z dv J> . 



t>b Jrn' ■ I \p 3 r 2 

The same argument, together with sharp versions of Povzner inequalities from 
E|, yields uniform bounds and appearance on every moments of the solution, as well 
as appearance of some exponential moments (this last point was first noticed in [2*5]). 
in a similar way as in [23 Proof of Proposition 3.2]. Indeed we prove the 

Proposition 3.1 Let g be a solution in C(R + ;L>1) H L 1 1 oc (IR + ; L3) to the rescaled 
Boltzman equation fli.i7| ) with e G [0, 1), with initial datum g in . Then it satisfies the 
following additional moment properties: 
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(%) For any s > 2, there is an explicit constant C s > 0, depending only on B, e, 
and g in , such that 

te[o,oo) 

(ii) Ifgm£ r ^ v G L 1 (R iV ) forr > andr] G (0, 1], there exists C\,r' > 0, depending 
only on B, e, and g in , such that 

sup f g{t,v)e r ' w dv <d. 
te[o,oo) Jr n 

(Hi) For any rj G (0,1/2) and t > 0, there exists a v ,C v G (0, oo), depending only 
on B , e, r and g- m , such that 

sup I g(t,v)e a ^ vlv dv <C V . 
te[T,oo)jRN 

Let us emphasize that the constant C s ,a^,C^ may depend on g m only through its 
mass p and its kinetic energy & m . 

Proof of Proposition \3.1l The proof is just a copy with minor modifications of 
classical proofs. For the proof of (i) we refer for intance to [HOI EU El] and the 
references therein. The proofs of (ii) and (iii) are variants of the proof of 
Proposition 3.2], which itself follows closely the proof of [5J Theorem 3] extended to 
the inelastic case in 9J. The starting point is the following differential equation on 
the moments 



d 

—m 
dt 



p = Qidi d) \v\ 2p dv + pm p with m p := / g\ 

JR N ' JR N 



\ 2p dv. 



Using the same notation as in [2B1 Proof of Proposition 3.2], we introduce the new 
rescaled moment function 



z. 



j " ~ T(ap+ 1/2)' ~ p ' fc=Wfc 



Z p :— max {zk+1/2 z P -k, z k z p-k+i/2}, 



for some fixed a > 2, and we obtain the differential inequality 



dy 

(3.2) ^ < A'p*/ 2 - 1 ' 2 Zp - A"p a / 2 z l p +1 l 2p +pz p 



for any p = 3/2, 2, . . . and for some constants A', A" > 0. Note that (|3.2|) is nothing 
but 28, equation (3.18)], with an additional term pz p due to the additional term 
— V„ ■ (v g) in equation ()1.17|) . 

On the one hand, we remark, by an induction argument, that taking p Q = 
Po(a, A', A") and x = x (a, A', A") large enough, the sequence of functions z p := x p 
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is a sequence of supersolution of (|3.2|) for any x > x and p > p . Let us emphasize 
here that we have to take a > 2 (i.e., rj < 1 in [2E1 Proof of Proposition 3.2]) because 
of the additional term pz p . On the other hand, choosing x\ large enough, which 
may depend on po, we have from (i) that the sequence of functions z p := x p is a 
sequence of supersolution of ()3.2|) for any x > X\ and for p G {0, 1/2, ...,p }. As a 
consequence, we have proved that there exists £2 := max(xo,£i) such that the set 

(3.3) C x := L = (z p ); z p < x p Wp G l - N 

is invariant under the flow generated by the Boltzmann equation for any x > x^'- if 
g(ti) G C x then #(£2) G C x for any £2 > £i- The end of the proof is exactly similar to 
that of |2HJ Proof of Proposition 3.2]. □ 

The integral upper bound in point (iii) of Theorem 11.31 follows from point (iii) 
of Proposition 13.11 

3.2 Stability in L 1 

The stability result [23 Proposition 3.3] translates for ()1.17j) into: 

\\g - h\\» + e~ 2T \\(g - h)\v\ 2 \\ L , < e c ^^ [\\g in - h in \\ L , + \\(g m - h m )\v\ 2 \\ L ,} 

for any solutions g and h in C(K. + , L\) nL°°(lR + , L3) with initial datum < g in , h in G 
L\. This shows that, in the Banach space L\, the evolution semi-group S t of ()1.17j) 
satisfies: for any t > 0, St is (strongly) continuous in any L\ bounded subset of 
L\. However we shall prove a more precise stability result, working directly on the 
rescaled equation (jl,17j) . 

Proposition 3.2 Let < g m , h in G L\ and let g and h be the two solutions of fll.l 7| ) 
(in C(M+, L\) fl L°°(1R + , L\) ) with e G [0, 1]. Then there is C sta b > depending only 
on B and sup 4>0 ||g + h\\ L i such that 



Vt > 0, \\g t - h t \\ L i < \\gm - h in \\ L i e 



stab 1 



Proof of Proposition ^. 21 We multiply the equation satisfied by (g — h) by 4>(t, v 
sgn(g(t, v ) — h(t, v)) (1 + \v | 2 ). We use on the one hand the same arguments as in j 
Proposition 3.4] to treat 

[Q(g,g) - Q(h,h)) <p(t,v)dv, 



which gives 

I<c([ (g + h) (1 + \v\ 3 )dv ) I / \g - h\ (1 + \v\' 2 ) dv 
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On the other hand we use that 



V„- (v(g-h))(f)(t,v)dv = -N \g - h\ (1 + \v | 2 ) dv 



+ / \g — h\V v ■ {v + v\v\ 2 ) dv 
2 \g — h\ \v\ 2 dv. 



This concludes the proof with C sta b = C sup t>0 \\g + h\\ L i +2. □ 

3.3 Uniform propagation of Lebesgue norms 

Let us take a normal restitution coefficient e G (0, 1) (the case e = can be included 
in dimension N = 3) and 1 < p < +oo, and let us consider some initial datum 
g- m G L\ fl L p . We compute the time derivative of the LP norm of the solution g to 
equation ()1.17|) : 

l*f gPdv= f Q + (g,g)g p - 1 dv- [ g?L(g)dv-[ g^ 1 V v {vg)dv. 

P at J R N J R N J R N J R N 

We use the control (jl.2Hj) . and 

/ V v -{vg)sr l dv = N(l--) \\g\\ p LP . 
Gathering all these estimates, we deduce 



~ / 9 p dv< [ Q + (g,g)g?- 1 dv 
P «r J r n j r n 



P 

Concerning the gain term, Theorem 12.21 yields, for any e > 



min<! 1,N ( 1- - 1 W g p (l + \v\) dv. 



[ Q + (g,9)f p - 1 dv<C £ \\g\\ 1 I t pe \\g 

Jm. n 2 



1+p9 """ P £- e) +e\\g\\ Li \\g(l + \v\)\\ p LP . 



Hence, using the bound Ce on the kinetic energy, if we fix e such that 

Cf-%<imin{l,iv(l-I 



< (1 llnll p(1 " 6 ' ) - K \\n\\ p 
\y\\LP — °+ \\9\\lp jx - \\9\\ l p 



we obtain 

d_ 

r iii/ -" Tii; " iw "-^ u ^ /P 

for some explicit constants C + , > 0. By maximum principle, it shows that the 
LP norm of g is uniformly bounded by 



sup \\gt\\ L P < max < [ -J- ] , ||#i n ||zp 



t>o 



The proof for weighted L p norms is exactly similar. This shows the part concerned 
with Lebesgue norms in point (i) of Theorem 11.31 
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3.4 Non-concentration in the rescaled variables and Haff 's 
law 

In this subsection we give a short proof of Haff 's law, even if a stronger pontwise 
estimate from below on the tail in rescaled variables will be proved in the next 
section. Let us take a normal restitution coefficient e G (0, 1) (the case e = can 
be included in dimension N — 3). Let f in = g in be an initial datum in L\ fl L p (with 
1 < p < +oo). Hence according to the previous subsection, the rescaled solution g 
to (jf .1 7j) with initial datum g in satisfies 

sup \\g t \\LP < C P 

t>o 

for some explicit constant C v > depending on the collision kernel and the mass, 
kinetic energy and L p norm of f- m . By using Cauchy-Schwarz inequality, this non- 
concentration estimate implies that for any r > 

Vt>0, I g(t,v)dv <Cr^ N . 

J\v\<r 

Thus there is r > such that 

Vt>0, / g(t,v)dv<% 

J\v\<r 

and thus 

(3.4) Vt>0, / g(t,v) \v\ 2 dv > 
Jr n 

> 
> 

As a conclusion, gathering (|3.f|) and (|3.4|) . we have proved that for some constants 
C , C\ G (0, oo) there holds 

Co<£{a{t r ))<c x , 

and Haff's law (|1.21|) follows thanks to (jl.2(J|) . which proves Theorem If .21 

Remark 3.3 The inequality £(f(t,-)) < M (1 + t)~ 2 (or equivalently £(g(t,-)) < 
C\) was already known: see for instance Q equations (2.5)- (2.6)] where it is proved 
for a quasi-elastic one- dimensional model with the same evolution equation hi. 9^ on 
the kinetic energy, by comparison to a differential equation. Indeed the harder part 
in Haff's law is the first inequality, which means that the solution does not cool down 
faster than the self-similar profile. As emphasized by the proof above, this is related 
to the impossiblity of asymptotic concentration in the rescaled equation fli.i7| ). 



g(t, v) \v | 2 dv 



\v\>ro 



r / g(t,v)dv 



\v\>ro 



g(t,v)dv > 



\v\<ro 



prl 
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3.5 Uniform propagation of Sobolev norms 

Let us take a normal restitution coefficient e G (0, 1) (the case e = can be included 
in dimension N = 3). The study of propagation of regularity and exponential 
decay of singularities is based on a Duhamel representation of the solution we shall 
introduce. Let us denote 



L{t,v) 



g{v*) \v — u*| dv* ) , 



and 



Stg = g(e t v) exp 
the evolution semi-group associated to 



■Nt- / L(s,e~ {t - s) v)ds 
Jo 



Tg 



Then the solution of (jl.l7|) represents as 



g(v*) \v - v*\ dv* ) g(v) - V v ■ [vg). 



gt = S t g in + / S t - s Q + (g s , g s ) ds. 
Jo 

We give a proposition similar to [HU Proposition 5.2]: 

Proposition 3.4 There are some constants a > 0, 5 > 0, K > and k > such 
that for any s, rj > 0, we have 



s t9m\\m +a ^ CWe Kt llflffcllff't? sup 



0<t<t 



s+k 



v+S 



St-sQ + (g s ,g s )ds 



< C Duh sup ||c/(V)||tf 



Jin 



0<t<t 



s+k 



v+S 



Proof of Proposition 13.41 The proof is exactly similar to 32, Proof of Proposi- 
tion 5.2]. Indeed the semi-group in [^21 Proof of Proposition 5.2] is 



St9 = g(v) exp 



Lis, v) ds 



and thus the estimates on the Sobolev norm in v can only improve for St according to 
St- The main tool of |H2l Proof of Proposition 5.2], i.e., the Bouchut-Desvillettes-Lu 
regularity result on Q + , has been proved in our case in Theorem 12.51 □ 



Now results follow as in 
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Theorem 3.5 Let < g m G L\ be an initial datum and let g be the unique solution 
of fll.l 7| ) in C{R + ,Ll) D L 1 {R + ,Ll) associated with g m . Then for all s > and 
T] > 1, there exists w(s) > (explicitly w(s) = 5\s/a], where a is defined in 
Proposition \3.4\ ) such that 



9m e if; 



ri+w 



m>\\g(t,-)\\m < +oo 

i>0 



with uniform bounds. 



Proof of Theorem \3.5l Let n G N be such that na > s (n = \s/a\). Let w(s) = 
8\s/ot\. The proof is made by an induction comprising n steps, proving successively 
that g is uniformly bounded in if ia ^ for i = 0, 1, . . . , n. 

Let us write the induction. The initialisation for i = 0, i.e., g uniformly bounded 
in L^ +w is proved by the previous study of uniform propagation of weighted LP norms 
in Subsection 13.31 Now let < i < n and suppose the induction assumption to be 
satisfied for all < j < i. Then proposition 13.41 implies 



< C*Duh e * 



|#m(-)llff«* . SUp ||#(W)|Dt-l)a 

V+^n, + S 0<t <t " V+ H=± w + S 



and 



S t -sQ + (gs,g s )ds 



< C Duh SUP ||^(* , -)ll^(t-l)a 

0<t <t , 1+ n=i w+5 



Moreover as i > 1, 



n — i n — (i — 1) 
r] H w + o < rj H w. 

n n 



Thus, using the induction assumption for i — 1, g is uniformly bounded in H %<x n _j , 
which concludes the proof. □ 



3.6 Exponential decay of singularities 

Let us take a normal restitution coefficient e £ (0, 1) (the case e = can be included 
in dimension N = 3). In this part we shall follow a similar strategy as in j22] in 
order to show that singularities decrease exponentially fast along the flow in rescaled 
variables. Namely we prove the 

Theorem 3.6 Let < g in G L \ R LP , p G (1, oo) ; and let g be the unique solution 
of in C(R + ,Ll) fl L 1 1 oc (lR + , L \) associated with g in . Let s > 0, q > be 

arbitrarily large. Then g can be written g s + g in such a way that 

sup ||^f |U |ni i < +oo, />0 

3A>0; \\gf\\ Ll =0{e- xt ). 

All the constants in this theorem can be computed in terms of the collision kernel, 
the mass and kinetic energy and L 2 norm of g- ui . 
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Proof of Theorem \3.b\ Assume first that < g- m G L\C\L P , p G [2, oo). Then g in G L 2 
and the proof of Theorem 13.61 is exactly similar to [321 Proof of Theorem 5.5] since 
the only tools of the proof are the stability result, the estimate on the Duhamel 
representation and the uniform propagation of Sobolev norms, which have been 
proved respectively in Proposition 13.21 Proposition 13.41 and Proposition 13.51 The 
propagation and appearance of moments in L 1 (used in this proof) were proved in 
Proposition 13. II Moreover as was already pointed out in (321 Section 7, Remark 3], 
it is possible with the same arguments to relax the assumptions on the initial datum 
to < g in G L\C\ L p for any p G (1, 2] by using the gain of integrability of the gain 
part of the collision operator. □ 

Point (i) of Theorem 11.31 is deduced from this theorem. 

Remark 3.7 We do not know how to carry the argument in VJ2\ Theorem 7.2] in 
order to reduce the assumptions to only < gi n G L\. The estimates of Abrahams- 
son on the iterated gain term can be easily extended to the inelastic framework, but 
the decomposition of Abrahamsson (quoted in the elastic case in I32\ Lemma 7.1]) 
between a part with finite L p norm for some p G (1, 3) and a part which decreases 
exponentially fast requires a lower bound on the energy. Here for the rescaled inelas- 
tic problem we deduce this lower bound from the propagation of L p bounds, which 
therefore seem compulsory in our method. 

Remark 3.8 A suggested by this study, the self-similar variables are not only useful 
for proving the existence of self- similar profiles, but it seems that they also provide the 
good framework for studying precisely the regularity of the solution. For instance, 
coming back to the original variables, Theorem 1 3. 61 shows the algebraic decay of 
singularities for the solutions of M.l)) . 

4 Self-similar solutions and tail behavior 

In this section we achieve the proofs of Theorem 11.11 and Theorem 11.31 by showing 
the existence of self-similar solutions, and obtaining estimates on their tail and the 
tail of generic solutions. We consider a normal restitution coefficient e G (0, 1) (and 
as before the case e = can be included in dimension N — 3). 

4.1 Existence of self-similar solutions 

The starting point is the following result, see for instance ITUl Theorem 5.2] or [3*lll8j. 

Theorem 4.1 Let y be a Banach space and (S t )t>o be a continuous semi-group on 
y. Assume that there exists JC a nonempty convex and weakly (sequentially) compact 
subset ofy which is invariant under the action of S t (that is S t y G fC for any y G K. 
and t > 0), and such that St is weakly (sequentially) continuous on JC for any t > 0. 
Then there exists y G fC which is stationary under the action of S t (that is S t yo = yo 
for any t > 0). 



30 



Proof of Theorem II . II (existence part). The existence of self-similar solutions follows 
from the application of this result to the evolution semi-group of (|1.17j) . The conti- 
nuity properties of the semi-group are proved by the study of the Cauchy problem, 
recalled in Sectional On the Banach space y = L\, thanks to the uniform bounds 
on the L\ and L p norms, the nonempty convex subset of y 

>c = l [ o<fey, f f( 1 v ) dv= (o) and II/IIlj + II/IIl,<m| 

is stable by the semi-group provided M is big enough. This set is weakly compact 
in y by Dunford-Pettis Theorem, and the continuity of St for all t > on K, follows 
from Proposition E21 This shows that there exists a nonnegative stationary solution 
to ()1.17|) in L\ fl L p for any given mass, that is a self-similar solution for the original 
problem (Jl.lj) . 

Then one can apply Theorem 13.61 which proves that the stationary solution 
of (jl,17|) obtained above belongs to C°° (in fact it proves that it belongs to the 
Schwartz space of C°° functions decreasing faster than any polynomial at infinity). 
Moreover, since the property of being radially symmetric is stable along the flow 
of ([1.170 . this sationary solution can be shown to exist within the set of radially 
symmetric functions by the same arguments. □ 



4.2 Tail of the self-similar profiles 

In this subsection we prove pointwise bounds on the tail behavior of the self-similar 
solutions. The starting point is the following result extracted from [§J Theorem 1]; 
notice that it is also a consequence of the construction of invariant sets C x for z p 
with a = 2, as defined in (j3.3[) . 

Theorem 4.2 (Bobylev-Gamba-Panferov) Let G be a steady state of fli.l7[ ] 
with finite moments of all orders. Then G has exponential tail of order 1, that 
is 

r* = sup < r > 0, / G(v) exp(r|t> |) dv < +oo > 
I Jr n ) 

belongs to (0, +oo). 

Note that if one define more generally (for s > 0) 

r* = sup < r > 0, / G(v) exp(r\v \ s ) dv < +oo 



a simple consequence of this result is that r* = +oo for any s < 1, and r* = for 
any s > 1. 

First let us prove the pointwise bound from above on the steady state. Since the 
evolution equation (jl,17[) makes all the moments appear (see Proposition 13. 1|) . we 
assume that G has finite moments of all orders. Moreover, as discussed above, we 
can also assume that G is smooth and radially symmetric. We denote r = \v\. We 
thus have the 



31 



Proposition 4.3 Let G G C 1 be a radially symmetric nonnegative steady state 
of \l.l r l\j with finite moments of all order. Then there exists Ai,A 2 > such that 

VveR N , G(v) < A ie ~ A2lvl . 

Proof of Proposition I4..':>l The differential equation satisfied by G = G(r) writes 

Q{G, G) — NG — rG' = 0. 

Since G is smooth and integrable, it goes to at infinity. By integrating this equation 
between r = R and r = +oo, we obtain 

G{R)=N rm dr -rmv ir . 

Jr r Jr r 

One deduces the following upper bound 

G(R)<N rm dT+ rQ-(^) dr . 

Jr r J R r 

Since Q~(G, G) = G (G * $), we have 

Q-(G,G)(v) <C(l + \v\)G. 
Hence, taking R > 1 leads to 

G(R) <C G(r) r N - x dr. 
Jr 

Finally, since we have by Theorem 14.21 

r+oo 

/ G(r) exp(A 2 r) r^" 1 dr < A < +oo 
Jo 

for some constants A ,A 2 > 0, we deduce that 

r+oo 

G{R)<C / G{r)r N - 1 dr<CA exp{-A 2 R) = A l exp{-A 2 R). 
Jr 

This concludes the proof. □ 

For the pointwise lower bound, we give here a proof based on a maximum princi- 
ple argument, inspired from the works [T^ I2(Jj. We shall in the next subsection give 
a more general result for generic solutions of (jl.l7|) , based on the spreading effect of 
the gain term and the dispersion (or transport) effect of the evolution semi-group 
of (jl.l7|) (due to the anti-drift term) in the spirit of [T3J|33]. 

Proposition 4.4 Let G G C 1 be a nonnegative steady state of {1.1 7| ) with finite 
moments of orders and 2 and which is not identically equal to 0. Then there exists 
Oi, a 2 > such that 

V«6R W , G(v) > aie-° 2H . 
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We first start with a lemma. 

Lemma 4.5 For any ro, a±, po, pi > 0, there exists ci2 > such that the function 
h(v) := ai exp(— a 2 \v\) satisfies 

(4.1) Vv,\v\>r , Q-(g,h)+V v -(vh) <0 

/or any function g such that 

g(v)dv = p 0l / gr(u) |u| = pi. 



Proof of Lemma \4. 51 On the one hand, it is straightforward that 

/i) := (o * $) h < (pi + po M) h. 
On the other hand, simple computations show that 

V v -(vh) = (N-a 2 \v\)h. 
Gathering these two inequalities there holds 

Vf, M>r , Q~{g,h)+V v -{vh) <(p 1 + N + p o \v\-a 2 \v\)h<0 

for a 2 large enough. □ 

Proof of Proposition \4.4\ Since G 6 C 1 and it is radially symmetric, there holds 
G'(0) = 0. As a consequence, the equation satisfied by G reads in v = 

Q(G, G)(0)-NG(0) = 

and then 

= Q+(G G)(0) 
K ' N 

since G is not zero everywhere. By continuity, G(v) > 2a± on 5(0, ro) for some 
ai,r > 0. 

Let us define 

Po := / G(v)dv, px : / G(t> ) |u| dv, 



and a 2 > associated to r ,di, p , pi by Lemma H31 On the one hand /i(t>) := 
ai exp(— a 2 \v\) satisfies (|4.1j) for g = G and, on the other hand, G satisfies 

(4.2) VveR N , Q~(G,G) + V v (vG) = Q + (G,G) > 0. 

Introducing the auxiliary function W := G — h, we deduce from (14.1)1 and ()4.2j) 

Vf , |u| > r , (G *$)W + V v (vW) > 

and VF(ro) = G(r ) — /i(r ) > G(r )/2 > 0. By the Gronwall Lemma (using that all 
the functions involved in this inequality are radially symmetric), we get W(v) > 
for any v, \v\ > Tq, which concludes the proof. □ 
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4.3 Positivity of the rescaled solution 

We start with three technical lemmas. 



Lemma 4.6 Let go satisfies forp G (1, oo) 

(4.3) [ g dv=l, [ g \v\ 2 dv<C 1 , / 9 P dv<C 2 . 

Jr n Jr n Jr n 

There exist R > r > and r\ > depending only on Ci, C 2 , and (vi)i=i,...,4 s^c/i i/ia£ 
\vi\ < R, i = 1, . . . , 4, and — fj| > 3r /or 1 < z ^ j '• < 3, and 



(4.4) 



/ g (v) dv>r] for i = 1, 2, 3, 

•/ B{vi,r) 



(4.5) Vtfj G B(vi,r), E e W3 WA n S 1 ^ W2 zs a sphere of radius larger than r, 

where E%, v stands for the plane defined in Proposition \1.5l and S^ Vt stands for the 
sphere of all possibles post-collisional velocity v' defined by 

Proof of Theorem \4.b\ Let Cr denotes the hypercube [—R,R] N centered at v — 
with length 2R > 0. Thanks to the mass condition and the energy bound in ()4.3|) . 
for R large enough, there holds 

(4.6) / g dv>± 

Then we define {Ki)%=\,...,i the familly of / = (2R/r) N hypercubes of length r > 
(with R/r G N), included in Cr and such that the union of Ki is almost equal to 
Cr. For any given A > to be later fixed, we may find r > such that 

(4.7) / g dv < \K t + 5(0, \r)\^' ( f g p dv) * < C [(A + l)r] N # < 1/4 

JKi+B(0,\r) \J Ki+B(0,\r) / 

for any % = 1, . . . , I. Hence we can choose K io such that the mass of go in Ki is 
maximal for % = i$. Because of ()4.6|) there holds 



(4.8) / g dv > 1/4(2 R/r 



-N 

y u,u ^ i/i^il/ / ) 

Gathering (|4.6|) and (|4.7|) we may find Kj Q C Cr such that dist(ifj , Kj ) > Ar and 
(14. 8|) also holds for i = j Q . 

Next, we fix A := 200/3. We define V\ (respectively v 2 ) as the center of the 
hypercube K io (respectively Kj ), and t>3 = (v\ + v 2 )/2 and V4 = v 2 . Then we have 

0(^3, ^4) = V 1 + — («2 - «l) e [fl, V 2 ), 
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which implies 

p-vi\ = ~y \v2~vi\ > -j- (Ar) > lOOr. 

Thus m fl S e (v 1 ,v 2 ) is a (N — 2)-dimensional sphere of radius larger than 100 r 
(because B(Q, lOOr) is included in the convex hull of S e (vi,V2)), and (|4.5jl follows 
straightforwardly. □ 

Lemma 4.7 Let us fix R > r > and rj > 0. Then there exists 5q > 0, rjo > 0, 
£o ^ (0, 1) (depending onR>r>0,r]>0 and B ) such that, for any functions f , h, 
£ satisfying \4-4\) - \4-^i) f or some velocities (ui)i=i,.„,4 such that \vi\ < R, i = 1, . . . , 4 
and \vi — Vj\ > 3r, 1 < i ^ j < 3, and for any £ G (£o, 1), i/iere holds 

Q + (f,Qt(h,£))>VolB(v 3 ,6 ), 
where we define here and below Q^(-, -)(v) = Q + (-, -)(£,v). 

Proof of Theorem \4.7\ We first establish a convenient formula to handle represen- 
tations of the iterated gain term. For any /, h and I and any v G M. N there holds 
(setting 'v = w and 'v* = u>*) 

Q + (f,Q+(h,£))(v) = C' b f l [ Q+(h,£)(w*)dw* \ dw. 

-V \ v - w \ [Jei w k J 

From the following identity 

Qf{h,£)(w*) = / / h(w 1 )£(w 2 )Qt{6 1 ,5 2 )(w*)dwidw 2 
Jr n Jr n 

where 5i stands for the Dirac measure at Wj, the term between brackets, that we 
denote by A, write 



A(v, w) 



j h{w l )£(w 2 ) <lim-3- / Q + (Si, <5 2 )(£ w*) H e (w*) dw* 1 dwi dw 2 

jR N xR N { £ ^ ^£jRN J 

where S e denotes the indicator function of the set {w* ; dist(«;*, E VjW ) < e}. Denot- 
ing now by D e the integral just after the limit sign in the term between brackets, 
and using the weak formulation ([1.250 . there holds 

£-1 r 

D £ = — — / 5i(z) 6i(z*) \z — z*\ b(a ■ z) S £ (£~ 1 z) dadzdz* 

^ £ JR N xR N xS N - 1 

= \w 1 -w 2 \t, C b da, 



N-l 



where in these integrates z' is defined from ( ) and next from {w\, w 2 , a) thanks 

to formula (jl.26j) . We define £o — (1 +r/R)~ l in such a way that z' — z'\ < r for 
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any z' G B(0,R) and £ G (£oj1)- Taking t> G B(v 3 ,r), w G Biv^r), W\ G £>(t>i,r) 
u; 2 G B{v 2) t), we have thanks to (|4.5j) and \w 2 — Wi\ > r: 



N-2 



D (v,w,w 1 ,w 2 ) := limD e >r£ 1 C b Cr 

£— >0 

As a consequence, for any t> G B(v 3 ,r), 

Q + (f,Qt(h,£))(v) > Q + (fl B{V4!r) ,Q + (hl B{vi>r) ,ei B{v2 , r) ))(v) 

> C' b [ [ [ — h(wx) £{w 2 ) D dw 1 dw 2 dw 

J B(v lt r)J B(v 2 ,r)J B(v 4 ,r) \V ~ W\ 

This concludes the proof. □ 
Lemma 4.8 For any v G M. N and 5 > 0, there exists k = k(5) > such that 
(4-9) Q + {v) :=Q + (1b(v,5)Ab(v,5)) >«l B ( flJ g,y 



Proof of Lemma 14.81 The homogeneity property (jl.l3|) of Q + and the invariance 
by translation allow to reduce the proof of ()4.9|) to the case v = and 5=1. The 
invariance by rotations implies that Q + is radially symmetric and the homogeneity 
property again allows to conclude that the support of Q + is a ball B' . More precisely, 
taking a C°° radially symmetric function such that > on I? = -8(0, 1) and < 
Is on R w , we have Q + (0, 0) is continuous, Q + > Q+(0, 0) on R N and Q+(0, 0) > 
on the ball B'. As a consequence, for any ball B" strictly included in B', there exists 
k > such that Q + > k1b«. In order to conclude, we just need to estimate the 
support of <2 + . 

Let us fix R G (0, 1) and choose 'v, 'v* G B(0, 1) such that 'v _L 'u*, | 'i>| = 
| = -R. Then for any a G S^ -1 , o _L 't> — the function Q + is positive at the 
post-collisional associated velocity v defined by 

'v + 'v* 1-e f , , 1 + e , 
v = — - — + — ^— ( v - v m ) + — ^— | v- v*\cr. 

Remarking that | 'v + 't>*| 2 — \'v — 'v*\ 2 = 2R 2 , ( 'v — 'v*) ■ {'v + 'f*) = and 
('v +'v*) • cr = \/2R, we easily compute 



v\ 2 = Rf 



1 + 



1 + 



>-R 2 , 
4 ' 



and the radius of B' is strictly larger than y/E/2. □ 
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Theorem 4.9 Let gi n satisfy the hypothesis of Theorem \l.cA and let g be the solution 
to the rescaled equation \1.1T\) associated to the initial datum g in . Then for any 
t* > 0, g(t, •) > a.e. on W N for any t > t*, and there exists ai,a,2,c> such that 

\/t>U, g(t,v) > aie~ a2|v| l M < ce *-** fora.e.veR N . 
Proof of Theorem \4.fA We split the proof into four steps. 

Step 1. The starting point is the evolution equation satisfied by g written in the 
form 

dtg + v- V v g + (N + \v\) g = Q + (g, g) + (\v\ - L(g)) g. 

Let us introduce the semi-group St associated to the operator v ■ V v + X(v), where 
\(v) := N + \v\. Thanks to the Duhamel formula and (jl.23j) . we have 

(4.10) g(t, •) > S t g(0, ■) + f S t ^ s Q + (g(s, -),g(s, ■)) ds, 

Jo 

where the semi-group S t is defined by 

(S t h)(v) = hive- 1 ) exp (- [ \{ve~ s )ds 

Notice that 



o 



- X(ve- S )ds^j > -(\v\ + Nt). 

Step 2. Let us fix t > and define g (t, •) := g(t + t,-). Using twice the Duhamel 
formula (|4.10|) . we find 



9o(t,-) > J S t ^ s Q + (g (s,.),J S s - s/ Q + (g (s',-),g (s',-))ds' J ds 

(4.11) > / / S t ^ s Q + (S s g ,S s ^ s ,Q + (S s >g , S s >g )) ds' ds. 

Jo Jo 

We apply now Lemma l4~o1 to g and set _R := 2R. Since St is continuous in L 1 , 
there exists Ti > 0, such that for any s G [0, Ti], there holds 

/ S s {g ){v)dv>r}/2 for 2 = 1,2,3, 

JB{v it r) 

and e~ Tl > £„. For v G 5(0, R ) and t G [0, Ti] we may estimate S t h from below in 
the following way 

(S t h)(v)> 1 h e - t (v) 

for some constant 7 = r fR 0t T 1 - The bound from below (|4.11|) then yields (using 
Lemma 14. 7|) 

Sto(V) > 7 2 J J Qt,^(s a go,Q^-S s »'9o,S 9 >goj S ) ds' ds 



> I 2 ri l ve s-t eB{v ^ r) ds' ds. 

Jo Jo 
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We have then proved that there exists T\ > and for any t\ G (0,71/2] there exists 
rji > such that (for some v G B(0,R)) 

Vt G [0, 71/2], •) := g (t + h, •) > Vl l B(Ml) . 

Step 3. Using again the Duhamel formula (|4.10|) and the preceding step we have 

9i{t,-)> / S t ^ s Q + (g 1 (s,-),g 1 (s,-))ds. 
Jo 

Thanks to Lemma f4. 71 on the ball -8(0, Rq), there holds 

gi{t, ■) > Vi S t - s Q + (l B (v,6i), 1b(v,Si)) ds 
Jo 

Jo 

> n 2 k($A p- ( - Ro+NTl h 1 /— 

^ Vl K {°l) e 1 1 B(i),v / T9<5i/4) 

on [0, T 2 ) with T 2 G (0, Ti/2] small enough, and then 

92(t,-) ■= 9i(t + > r) 2 l B (v,s 2 ) on [0,T 2 /2] 

with <5 2 := a/I9<5i/4 and t 2 G (0,T 2 /2] arbitrarily small, r\ 2 > 0. Repeating the 
argument we obtain 

9k{t, == ^ + Xj tiJ ") - % 1b ^) on t°' T ^/ 2 ]' with 5 * := (v / 19/4) fc *i 

with k > 1 and some tj G [0, Tj/2] arbitrarily small, rjk > 0. As a consequence, 
taking k large enough in such a way that <5fci?o, we get for some explicit constant 
1]* > and some (arbitrarily small) time t* > 



(4.12) Vt >0, <?(t* + t ,-) >^1 



B(0,R)- 



Step 4. Coming back to the Duhamel formula ()4.10|) where we only keep the first 
term, we have, for any to > 0, 

Vt > U, g(t + t,v) > r)*l\ v \< Re t-u exp(-|u| -N(t-Q). 

As a consequence, for any t > t*, 

(4.13) g(t,v) > l\ v \< Re i-t* [ sup l\ v \ =Re s exp(-|t>| - N s) 

\se[o,t-u] 

> l\v\<Re*-t* sup l\ v \ =Re s exp(-|u| -JV ln + (\v\/R)), 

\se[0,t-t*] J 

and we conclude gathering ()4.12|) and ()4.13|) . □ 

It is straightforward that Theorem 14.91 implies the lower bound in point (iii) of 
Theorem 11.31 
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5 Perspectives 



As a conclusion, we discuss some possible perspectives arising from our study (partial 
answers to them shall be studied in a forthcoming work [29 ). 

Let us denote 

V = |g G C°°, G radially symmetric, 

3a 1 ,a 2 ,A 1 ,A 2 >0 | a ie " a2H < G(v) < Ae"^ 1 }. 

Conjecture 1. For any mass p > 0, the self-similar profile G p with mass p and 
momentum is unique. 

If Conjecture 1 is true, the natural conjecture is 
Conjecture 2. (Strong version) For any initial datum with mass p and momen- 
tum (maybe with some regularity and/or moment assumptions), the associated 
solution satisfies (in rescaled variables) 

9t ~*t-*oc G p , 

where G p is the steady state of (|1.17j) with mass p and momentum 0. 

A relaxed version can be 
Conjecture 2. (Weak version) For any initial datum with mass p and momentum 
the associated solution satisfies (in rescaled variables) 

g t = g? + g? 

with gf EV and gf 1 -^t~-><x in L 1 . 

Note that the weak version of Conjecture 2 still makes sense when the self-similar 
profile with mass p and momentum is not unique and even if there is no convergence 
towards some self-similar profile (which could be the case for instance if the solution 
in rescaled variables "oscillates" asymptotically between several self-similar profiles). 
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